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Abstract. Delay-tolerant networks (DTNs) are characterized by a possible absence of end-to-end communica- 
tion routes at any instant. Still, connectivity can generally be established over time and space. The optimality 
of a temporal path (journey) in this context can be defined in several terms, including topological (e.g. shortest 
in hops) and temporal (e.g. fastest, foremost). The combinatorial problem of computing shortest, foremost, and 
fastest journeys given full knowledge of the network schedule was addressed a decade ago (Bui-Xuan et al, 
2003). A recent line of research has focused on the distributed version of this problem, where foremost, shortest 
or fastest broadcast are performed without knowing the schedule beforehand. In this paper we show how to build 
fastest broadcast trees (i.e., trees that minimize the global duration of the broadcast, however late the departure 
is) in Time- Varying Graphs where intermittent edges are available periodically (it is known that the problem is 
infeasible in the general case even if various parameters of the graph are know). We address the general case 
where contacts between nodes can have arbitrary durations and thus fastest routes may consist of a mixture of 
continuous and discontinuous segments (a more complex scenario than when contacts are punctual and thus 
routes are only discontinuous). Using the abstraction of T-CLOCKS to compute the temporal distances, we solve 
the fastest broadcast problem by first learning, at the emitter, what is its time of minimum temporal eccentricity 
(i.e. the fastest time to reach all the other nodes), and second by building a foremost broadcast tree relative to 
this particular emission date. 



1 Introduction 

Highly-dynamic networks, in particular Delay-Tolerant Networks, are characterized by a possible absence of end- 
to-end communication routes (paths or direct journeys) at any time. In most cases, however, communication can 
still be achieved over time and space using store-carry-forward-like mechanisms over disconnected routes (indirect 
journeys). This singularity led researchers to develop a range of routing techniques based for example on prior 
knowledge of the dynamics B4I16I . probabilistic information [23 1, encounter-based strategies 0101131171 . or contact 
statistics |20|. On the other hand, considering the time-dimension led to the extension of many graph theoretic and 
analytical concepts including paths and reachability H2I181 , distance |4 15 1, diameter [9 |, tree width (25], connec- 
tivity H 1131 , or necessary conditions [5. 22]. Of particular interest in this paper are the concepts of journeys (temporal 
paths) and temporal distance. 

The fact that connectivity takes place over time implies that the optimality of a given journey does not depend 
anymore on the sole number of hops separating the nodes. Short journeys might have a long duration, and long 
ones be comparably fast. In fact, at least three optimality metrics could be considered for a journey, as pointed 
out in |4): being shortest (least number of hops), foremost (earliest arrival date), or fastest (minimum time spent 
between departure and arrival, however late the departure is). The problem of computing shortest, foremost, and 
fastest journeys was addressed in [4| in the context of centralized algorithms using full prior knowledge of the 
network schedule. 

In a recent line of research 116171 , the authors started to look at the distributed variant of this problem, namely 
performing foremost, shortest, and fastest broadcast in time-varying graphs without knowledge of the schedule. 
Specifically, we asked what assumptions could make each of these broadcasts possible, and what others would 
make the broadcast trees reusable for subsequent broadcasts. Within the set of assumptions considered, the easiest 
problem regarding feasibility turned out to be foremost broadcast, while the easiest regarding reusability is shortest 
broadcast J6). Due to its intrinsic features, fastest broadcast remained unfeasible in all the contexts considered, 



which motivates us to consider stronger assumptions such as periodicity. The periodic assumption holds in practical 
scenarios such as transportations e.g., where entities have periodic movements (satellites, trains, or buses). Previous 
works in periodic DTNs include exploration 111 11121141 or scalable routing Ml 91241 . 

We are tackling the problem in two steps. First, having the emitter decide when broadcast has the potential to 
be the fastest, then build a foremost broadcast tree relative to (any of) the corresponding emission date. We show 
that the latter is trivial, and thus mainly focus on the way the emitter could process the optimal date to initiate 
a broadcast - i.e., its moment of minimum temporal eccentricity. This processing is achieved by using a recently 
introduced primitive, T-CLOCKS, which allows nodes in an arbitrary dynamic network to learn temporal lags a 
posteriori (that is, from a receptor-based perspective). Temporal lags relative to the emitter are thus first learnt by 
destination nodes {i.e., all the other nodes), then aggregated back to the emitter in an opportunistic way so that the 
emitter learns the complete evolution of its temporal eccentricity over one period p. Since the network is periodic, 
this information suffices to determine the emitter's temporal eccentricity at any time in the future. One of the main 
assets of our solution is that it addresses the general (and more realistic) case where contacts between nodes can 
have arbitrarily durations and possibly overlap in time with other contacts. This feature renders temporal distances 
substantially more complex to process, because they involve the co-existence of continuous and discontinuous 
routes in the network (a vast majority of DTN-related results assume punctual contacts and thus can only deal with 
discontinuous routes). 

The paper is organized as follows: in Section [2] we describe the model and terminology used in the paper. 
Then Section |3]present the building block upon which our algorithm relies, that is, the temporal-lags vector clocks 
(T-CLOCKS). Finally, Section [4] describes the algorithm, proves its correctness, and illustrate it with an example 
scenario. 

2 Background 

2.1 Model and assumptions 

Consider a set V of n nodes, making contacts with each other over a (possibly infinite) time interval T C T, called 
lifetime of the system; the temporal domain T corresponds here to M + (continuous-time system). Let the contacts 
between nodes define a set of intermittently available undirected edges E C V 2 such that (x,y) e£tti and y 
interact at least once in T. 

Following [8 1, we represent the network as a Time-Varying Graph (TVG, for short) Q — (V, E, T, p, C), where 
p : E x T — > {0, 1}, called presence function (or schedule), indicates whether a given edge is available at a 
given time, and ( e T, called latency, indicates the time it takes to propagate a message over an edge. In this 
work, we assume the latency is constant for all edges and presence times (though in general it could be defined 
as a function) and known to the nodes. For this paper, we consider Periodically-Varying Graphs (PVG, for short) 
where intermittent edges are available periodically. That is, the schedule p is unknown to the nodes, but assumed 
periodic: Ve G E, Vt G T, Vfc G N, p(e, t) = p(e, t + kp), and the period p is known to the nodes. If a message is 
sent less than C, time units before the disappearance of an edge, it is lost. Given an edge e, we allow the notation 
P[ti,t 2 )( e ) = 1 to signify that Vt G [ti, tq), p(e, t) = 1 (i.e., the edge is continuously present during interval [t\, £2))- 

Given a time-varying graph Q — (V, E, T, p, C), we denote by G = (V, E) its underlying graph, that is, in a 
sense, the static counterpart of Q. A sequence of couples J — ((ei, ii), (e 2 , £2) • ■ • , [&hi tk))< where e±, e2, 
is a walk in G, ti + £ < i J+ i and P[ti,ti+()( e i) = 1 f° r all 1 < i < fc is called a journey in Q. We denote 
by departure{J), and arrival(J), the starting date t\ and last date tk + C of J, respectively. Journeys can be 
thought of as paths over time from a source to a destination and thus have both a topological and a temporal lengths. 
The topological length of J is the number \ J\h = k of couples in J (i.e., number of hops), and its temporal length 
is its duration \J\ t — arrival(J) — departure(J) = tk — t\ + (,. For example the journey ((ac, 2), (cd, 5)) in 
Figure[T|has a topological length of 2, and a temporal length of 3 + £ time units. 

Let us denote by Jg the set of all journeys in time-varying graph Q, and by J* u ^ G Jg those journeys starting 
at node u and ending at node v. 

We say that a journey is direct if the presence of every two successive edges overlap in time and their use 
follow on directly; it is said indirect otherwise. An example of direct journey in the graph of Figure Q] is J\ — 
{(ab, 2), (be, 2 + ()}. Examples of indirect ones include J 2 = {(ac, 2), (cd, 5)}, and J 3 = {(ab, 2), (be, 2 + 

C),M,5)}. 
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Fig. 1. A time-varying graph Q. {Labels indicate when an edge is present.) 



We assume nodes can instantly detect the appearance or disappearance of an incident edge, and associate ded- 
icated operations in reaction to these events. Since the presence intervals are by convention right-open, we assume 
that the disappearance of an edge at time t is always handled before the appearance of another edge at date t, 
which is consistent with looking at journeys whose presence intervals strictly follow each other as indirect ones 
(e.g. J( a ,d) — {(ac, 5— £), (cd, 5)} in the above example). Processing times are assumed negligible. The entities 
need not sharing a common global time, but their clocks must advance at a same rate. Finally, we assume nodes 
have unique identifiers. 

2.2 Definition of the probiem 

As mentioned above, the length of a journey can be measured both in terms of hops or time. This gives rise to two 
distinct definitions of distance in a graph Q: 

- The topological distance from a node u to a node v at time t, noted d u j(v), is defined as Min{\J\h : J G 
J* u . A departure(J) > t}. For a given date t, a journey whose departure is t' > t and topological length is 
equal to d u ,t{v) is called shortest ; 

- The temporal distance from u to v at time t, noted d u .t{v) is defined as Min{arrival(J) : J € J* uv ) A 
departure(J) > i} — t. Given a date t, a journey whose departure is t' > t and arrival is t + d u ,t{v) is called 
foremost, and one whose departure is t' > t and temporal length is Min{d u ,t"{v) : t" > t} is called fastest . 

Informally, a fastest journey is one that minimizes the time spent between departure and arrival (however late the 
departure is). This metric is particularly relevant in communication networks whose medium is shared exclusively 
(e.g. to minimize the holding time), or transportation networks (to minimize a trip duration, even if this implies leav- 
ing later). The problem of computing shortest, fastest, and foremost journeys in delay-tolerant networks was solved 
in @ as a centralized (i.e., combinatorics) problem, given complete knowledge of Q. We consider a distributed 
variant of the problem, namely that of performing fastest broadcast, which we define as follows. 

Fastest Broadcast: Given an emitter src and a message m, m should be sent from src to all the nodes in the 
network in such as way that the global duration between the first message emission (at src) and the last message 
reception (anywhere) is minimized. We also require that the emitter detects the termination, whether implicitly or 
explicitly, however this detection does not need to be the fastest. 

2.3 Working with temporaf views 

A central concept in this paper is that of temporal view, introduced in [21 1 in a context of social network analysis. 
(This concept was simply called "view" in ||2TI . but we added the "temporal" adjective to avoid confusion since 
"view" has a different meaning in distributed computing, e.g. l27l .) The temporal view a node v has of a node u at 
time t, denoted <t> v ,t{u), is the latest (i.e., largest) t' < t at which a message received by time t at v could have been 
emitted at u; that is, in our formalism, 



There is a clear connexion between temporal distances and temporal views. Both actually refer to the same 
quantity seen from different perspectives: the temporal distance is a duration defined locally to an emitter at an 
emission date, while the temporal view is a date defined locally to a receptor at a reception date. In fact, we have 



4>v,t(u) = Max{departure(J) : J e JL V ) A arrival(J) < t}. 



d u ,te (v) =t r - (j) v j r (u) 



(1) 



where t e is an emission date, and t r is the corresponding earliest reception date. 

This definition can be seen as a generalization of the one in ||2TI in which only punctual contacts were con- 
sidered. As pointed out in |6|, the temporal view a node has of another is deeply impacted by the co-existence of 
direct and indirect journeys. Indeed, assuming arbitrary long contacts between nodes makes it possible for adjacent 
contacts to overlap in time and thus produce more complex patterns of time lag between nodes. Consider the plots 
in Figure |2] showing an example of evolution of a temporal distance (from a to c) and the corresponding temporal 
view (that c has of a) in a very simple TVG. Contrary to the case with punctual contacts - where evolution occurs 
only in discrete steps - there is here a mixture of discrete and continuous evolution. (The reader is encouraged to 
spend a few minutes on this example as these concepts are essential in what follows.) 
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(a) A simple TVG (b) Temporal distance from a to c (c) Temporal view that c has of a 
Fig. 2. Temporal distance and temporal views as a function of time (with £ 1). 



Direct journeys are often faster than indirect ones, but this is not necessarily the case (imagine a very long direct 
journey, versus a short indirect one whose edges traversals follow closely). As a result, the temporal view at anytime 
could be caused by either types of journey. Let us call direct view and indirect view the views resulting from the best 
direct and indirect journeys at any given time (keeping in mind that the temporal view is always the max between 
both). 

Direct and indirect views are of a different nature. Indeed, contrary to indirect journeys, which take place in a 
discrete way, direct journeys in general belong to a continuum of several such journeys, which induce continuous 
increases of the view. Looking again at the example on Figure [2(c)| the (direct) view c has of a during [1 + £, 3) only 
depends on the topological length of the corresponding journey, here 2 (node c can receive any message emitted by 
a between times 1 — £ and 3 — 2£ after a lag of exactly 2( time). Since the duration of a direct journey depends only 
on its topological length, it is sufficient for a node u to know exactly what is the length of the shortest direct journey 
currently arriving to it from another node v, called the level of u with respect to v, to deduce the corresponding 
direct view in real time. This particular way of defining a level can be formalized as follows: 

level Vlt (u) = Min{\J\ h : J £ JL V \ A isDirect(J) A arrival(J) = t}, 

where isDirect(J) is true iff J is a direct journey. By convention, level v> t(u) is considered to be — oo if no 
direct journey from u is arriving to v at time t. Let us insist on the fact that our notion of level is relative to the 
reception side, and is not concerned for example with the fact that some edges of the journey may have disappeared 
by reception time. What matters is to know whether messages could still be currently arriving from a given remote 
node through a given number of hops. 



3 Tracking temporal views with T- Clocks 



This section is concerned with tracking in real-time the evolution of direct and indirect views. A powerful ab- 
straction in this regard was provided in [7 |. Precisely, an algorithm called T-CLOCKS was provided that allows 
to continuously keep track of direct and indirect views in a delay-tolerant network, and allows to plug another 
algorithm on top of it to solve more concrete problems using this information. The example given in |7 j was the 
construction of foremost broadcast trees in periodically-varying graphs. In this paper, we will use T-CLOCKS to 



solve the construction of fastest broadcast trees in periodically-varying graphs, which uses a completely different 
approach compared to the foremost case. 

Concretely, T-CLOCKS maintains two kinds of variables: the level of the local node with respect to any other 
node (as previously defined, which accounts for the corresponding direct view), and the largest date at which a 
message carried to the local node through an indirect journey could have been emitted at that other node (accounting 
for indirect views). Note that even though the nodes clocks might not be synchronized, the fact that they all advance 
at the same rate and the latency is known to the nodes allows T-CLOCKS to deliver information in the local time 
referential. 

Technically, the abstraction materializes as an intermediate layer between the network and some higher algo- 
rithm (see Figure [3}, which is informed by means of generating the two following events: levelChanged (src, 
level, proxy), reflecting the evolution of a direct view, and date Improved ( src, date, proxy ) , reflecting 
that of an indirect view. (Both events are further described below.) 




Network 



Fig. 3. T-CLOCKS as an abstraction to track temporal views. 



Given a direct or indirect view, we call proxy the local neighbor responsible for making the view evolve, i.e., 
the last node before the local node in the corresponding journey (this information is crucial to learn routing paths 
in the network). Notifications do not only occur when a level or a date has been updated, but also occur when a 
proxy has changed (which may happen without any change of view, e.g. when two local neighbors are providing the 
same direct view relative to a given remote node, and one of them disappears). T-CLOCKS are used in the spirit of 
the observer-observable design pattern for Object-Oriented development [26] between two objects, namely the T- 
Clocks algorithm (observable) and the higher algorithm (observer). The higher algorithm subscribes to the events 
of the T-CLOCKS algorithm by calling a registerQ function. Then it is notified of both types of events by means of 
calls to two functions levelChangedQ and datelmprovedi) to be implemented on the higher algorithm and that 
encapsulate the desired response to the events. The notifications are raised as follows: 

- levelChanged(Node src, Integer level, Node proxy): called whenever the proxy or direct view with respect 
to another node (also called source) has changed. This information is thus received as parameters of the call. 

- dateImproved(Node src, Integer date, Node proxy), called whenever the indirect view relative to another 
node (src) has increased if and only if this date is larger than the direct view (that is, larger than now() — 
level x £ for the corresponding source). 

T-CLOCKS are typically assumed to run independently from the higher algorithm, and we will assume that it 
is already running when our higher algorithm registers to it. Note that using T-CLOCKS does hide a substantial 
amount of complexity to the higher algorithm, and indeed the additive cost of our fastest broadcast tree construc- 
tion algorithm is marginal. The complexity of T-CLOCKS is unknown as of today, but the motivation is great to 
characterize and improve it, as any such improvement would have direct repercussions on the performance of all 
algorithms that use it as building block. 



4 Learning fastest broadcast trees in periodic TVGs 



We now describe how to build fastest broadcast trees in periodically-varying graphs. This problem can be tackled in 
two steps. First, having the emitter learn when the broadcast has the potential to be the fastest, then build a foremost 
broadcast tree relative to any of the corresponding emission dates (by definition, nothing better can be done than 
a foremost broadcast tree once the best dates are known). Building a. foremost broadcast tree relative to a given 
emission date does not require specific or difficult processes. This can be done by means of a flooding where all 
the nodes record which of their neighbor gave them the message first, followed by local acknowledgments of these 
relations. Thus, the rest of the section focuses on the mechanisms by which a node can learn when the broadcast 
has the potential to be the fastest {i.e., its moment of minimum temporal eccentricity) using T-Clocks. 
The temporal eccentricity (or simply eccentricity below) of a node u at date t is formally defined as 

ecc u (t) = max{d u , t (v) : v G V}, (2) 

that is, the maximum among all temporal distances (or simply distances below) from u to all other nodes at time t, 
see Section |Z21 for definitions. Informally speaking, it is the duration of the "slowest" foremost journey departing 
at t from node u. 

4.1 High-level informal strategy 

The algorithm consists in inferring (and recording) temporal distances at every node relative to a given emitter, 
based on the evolution of temporal views monitored through T-Clocks. Precisely, for a given emitter u, every 
node v ^ u infers d u (v) from 4> v (u) over one period and records it in a table called distance table. Since we 
deal with continuous-time and possibly overlapping contacts, this information is recorded as a set of intervals that 
correspond to the different phases of evolution of the distance (discrete or continuous). The distance tables of all 
nodes are then opportunistically aggregated (arbitrarily) along a tree rooted in u. The aggregation of a children table 
consists of a segment-wise maximum against the local distance table (segments can be artificially split to match). 
Once the emitter has aggregated the table of its last child, the final result corresponds to its eccentricity over time 
(Equation |2J. It finally selects any of the minimum values as preferred initiation date for the intended broadcast, 
then terminates. 

4.2 Learning the temporal eccentricities over a period 

Let us first observe that learning the minimum temporal eccentricities in distributed networks is an interesting 
problem in its own right, and even though we use it here as a primitive for broadcasting, it could be used as well for 
other tasks, such as electing a leader based on its ability to reach all others quickly. 

We now describe how tables of distances are learnt on the receptor side using T-Clocks. Based on the re- 
lationship between temporal view and temporal distance, we establish three key properties through the following 
Lemmas. In what follows we use interchangeably the terms emission date and initiation date, both referring to the 
time when a potential message is sent at the emitter. 

Lemma 1. Every discrete increase of(j> v (u) by value k corresponds to a continuous decrease of d u {v) of duration 
k. 

Proof. A discrete increase by k of the view at time t (that is, (j) v .t(u) = </>„. t _ e (u) + k), implies the existence of 
a journey J from u whose departure is and arrival is t. Switching to the emitter viewpoint, the consid- 

ered increase implies that no journey starting during [ti = (f> Vi t-e(u), t<i = <fi v .t(u)) could have arrived before t 
(otherwise such a journey would imply an intermediate increase of the view by less than k). Therefore, the tem- 
poral distance at any point in [ij., is) is fully determined by the arrival of J, and thus for all t' in [t\, ti), we have 
d u ,t' (v) = d u ,t 2 ( v ) + (*2 — t'), which corresponds to a continuous decrease of the distance during t^ — ti^k time 
units. □ 

Lemma 2. Each continuous increase of <f> v (u) during k time units corresponds to a stagnation of d u {v) during k 
time units. 



Proof. Continuous increases of the view are due to continuums of direct journeys of same level. Such increase 
during some interval [t, t + k) thus implies a continuum of direct journeys departing over [(j> Vtt (u), <f> Vt t(u) + k) 
(since £ is constant). We thus have Vt e e [<Au t(u), 4>v.t( u ) + k), d u j e (v) — level x £ (where level is the level of 
the considered journeys), which corresponds to a constant during k time units. □ 

Lemma 3. All initiation dates are covered either by a discrete or a continuous increase of the view. 

Proof. By nature of the view which is past-inclusive, i.e., the fact that a message emitted at time t e could have 
arrived by some time t r implies that any message emitted before t e could have also arrived by t r . (In essence, there 
is no "gap".) Besides, an increase is necessarily discrete or continuous. □ 

Combination of Lemmas[T][2]and[3]allows us to state the following general theorem on temporal distances (with 
constant edge latency Q. 

Theorem 1. The evolution of d u (v) can be fully captured by a sequence of segments of two possible types: flat 
segments (stagnation of the value) and slope segments (continuous decrease of the value). This sequence can be 
inferred at v by associating every continuous (resp. discrete) increase of <j) v (u) to a flat (resp. slope) segment of 
d u (v). 

This strategy is the one considered by AlgorithmQ] whose underlying principle consists in detecting (and record- 
ing) every transition between two segments of distance by means of transitions in the evolution of the temporal view 
(using T-Clocks). Each transition is recorded as a triplet containing i) an emission date, ii) the corresponding value 
(distance), and iii) the type of segment starting at this emission date (flat or slope). An example of such sequence 
is given on Figure[5] representing the distance from node a to node c in the example TVG of Figure|4] 
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Fig. 4. Example of periodic TVG (with period p = 100 and edge latency ( — 1). 
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(a) Sequence of triplets for d a ,t (c) „ f a ta m ta — 1ST 

(b) d a ,t(c) as a function of t 

Fig. 5. Temporal distance from a to c, as a function of the emission date. 



The subtlety here is that the value of the distance at the beginning of a slope segment, as well as the duration 
of both types of segments, becomes known only at the end of the corresponding temporal view segment; therefore, 
the algorithm always records values relative to a previously pending emission date (pendingED). Precisely, when- 
ever an event related to the temporal view occurs, whether it be caused by the arrival of a better indirect journey 
(datelmprovedQ) or a change in the level (levelChangedQ), the same update function is called involving the 
following sub-routines (see Algorithm[T): 



Algorithm 1 Computing temporal distances at a receptor, relative to a given emitter. 



1: Map <Date, Distance> table <s— 

2: Integer currentLevel < hoo 

3: Date startD <— nil 

4: Date pendingED <— nil 

5: initQ: 

6: TC locks. register^) 

7: datelmproved (Node src, Date date, Node proxy): 
8: if src = emitter then 
9: update(date) 

10: levelChanged(Node src, Level level, Node proxy): 

11: if src = emitter then 

12: update(now() — level x 

13: currentLevel «— ieuei 

14: update (Date newED): 

15: it startD — nil then 

16: startD <— now() 

17: pendingED <— newED 

18: else 

19: updateFlatQ 

20: updateSlope(newED) 

21: if now() = startD + p then 

22: terminate 

23: updateFIat (): 

24: if currentLevel < +oo then 

25: Date bestED <— noiu() — currentLevel x £ 

26: if bestED > pendingED then 

27: D> Aflat segment is detected 

28: table.add(pendingED, currentLevel x C,," flat") 

29: pendingED <- bestED 

30: updateSlope (Date new ED): 
31: it new ED > pendingED then 
32: D> A .s/ope segment is detected 

33: table.add(pendingED, now () — pendingED, "slope") 
34: pendingED <— newED 



- updateFlatQ: If the current level (currentLevel) was not infinite, then the corresponding continuum may have 
delivered new emission dates (checked in lines l24ltol2oTi. If this is the case, then a flat segment is inserted for the 
pending date using distance value currentLevel x £ (see proof|2]i. 

- updateSlopeQ: If the new event implies a discrete improvement of the view (line l3TTi. then a slope segment 
must be created, and only then the distance at the pending emission date becomes known. The segment being 
a continuous decrease, the value corresponds to the distance at the newly received emission date plus the time 
elapsed between that date and the pending date (see proof[T]l, i.e., (now() — newED) + (newED — pendingED) = 
now() — pendingED (see line[33]l. 

Observe that nothing prevents a same event from inducing both a flat and a slope segments, as is the case in the 
triangle TVG of Figure [4] when c's level relative to a changes at time 21. Table Q] shows an example of execution 
trace corresponding to the most relevant steps at node c (relative to emitter a) in the example of Figure [5] The 
beginning of execution at node c was arbitrarily set anytime between date 21 and date 60 (modulo p = 100). 



Date 


Event 


Action 


60 


levelChanged(+oo) 


startD <- 60 
pendingED <s— 59 
currentLevel < — hoo 


111 


lev elC 'hang ed(2) 


updateFlatQ 
updateSlope(109) 
table.add(59, 52, "slope") 

pendingED <— 109 
currentLevel <— 2 


121 


lev elC hang ed(l) 


updateFlatQ 

table.add(109, 2, "flat") 

pendingED ^ — 119 

currentLevel <— 1 

updateSlope(12Q) 

table. add(HQ, 2, "slope") 

pendingED <— 120 
currentLevel <— 1 


160 


levelChanged(+oo) 


updateFlatQ 

table. add(120, 1, "flat") 

terminate 



Table 1. Relevant traces at node c relative to emitter a (assuming an arbitrary start of execution between dates 21 and 60). 



Theorem 2. The execution of Algorithm Q] in a periodically-varying graph Q with known period p, relative to 
emitter u, results in every node v capturing correctly the evolution of d u (v) over one complete period, in 0{p). 

Proof. The argument is based on the correctness of T-CLOCKS [7| combined with the strategy of TheoremQ] which 
Algorithm Q] implements. Theorem[T]states that each segment of evolution of the temporal distance corresponds to 
a segment of evolution of the temporal view (the precise characterization of this correspondence being given by 
the proofs of Lemma Q] and |2). The correctness follows from detecting all transitions in evolution of the temporal 
view (events levelChangedQ and datelmprovedQ), guaranteed by T-CLOCKS, and transposing the correspond- 
ing segments into segments of the evolution of the temporal distance by means of procedures updateFlatQ and 
updateSlopeQ in Algorithm Q] (both of which are called for every such transition, and each of which checks for 
the need to create the corresponding type of segment according to Lemmas Q] and |2j. One period exactly after the 
first record is inserted, a last record is inserted in the distance table, then the algorithm terminates (linesl2Tlandl22l. 
Because the dates are taken modulo p, this last record completes the distance table relative to one complete period. 

□ 

4.3 Aggregating distance tables back to the emitter 

Distance tables relative to a given emitter are aggregated along a tree rooted at that node. Such a tree may be arbitrary 
and built, for instance, using the same strategy as mentioned above for single foremost broadcast trees (i.e., flooding 
a dedicated message from the emitter and detecting from which neighbor this message is first received at every node, 
followed by local acknowledgments of the corresponding relations). Once a node has computed its distance table 
and aggregated the tables of all its children (if any), it sends the resulting table to its parent. The aggregation of 
a children table (described below) consists of a segment-wise maximum among both tables for all emission dates 
(see Figure |6]for a visual illustration). Once the emitter has aggregated all its children tables within its own table 
(initially consisting of a single flat segment of value 0), the final result corresponds to the maximum distance in 
every point in time among all nodes, and therefore (Equation [2]) represents the evolution its eccentricity over one 
period. Once this table known, it finally selects any of the dates at which the eccentricity is minimum (relatively to 
which the foremost broadcast tree is to be built), then terminates. We now describe the aggregation process in more 
detail. 

The purpose of aggregation is to select the maximum value among all distances in every point in time. Aggregat- 
ing two tables based on their segments would be relatively straightforward if the segments were horizontally aligned 




(a) d a ,t(c) as a function of t (b) d a%t (b) as a function of t (c) ecc a (t)=max(d a ,t(b) , d a ,t(c)) 

Fig. 6. Aggregation of distances. The left curve (distance from a to c) is combined to the middle curve (distance from a to b) in 
order to yield the eccentricity of node a over one period. 

(i.e., if the period was split into the same sequence of intervals in both tables) and vertically non-crossing (i.e., given 
any pair of aligned segments, one of them remains higher than or equal to the other during the corresponding inter- 
val). In such an ideal case, the maximum operation between two aligned segments, say <Sj = (U,d ti , trendi) and 
S'i = (ti, d' t . , trend'i) is straightforward; it may consist of selecting <Sj iff 

(d ti > dt'.) V (d ti = dt 1 . A trendi = flat); and S[ otherwise. 

Theorem 3. Given two segments whose intervals are aligned and values are non-crossing, the above logic selects 
the correct maximum segment. 

Proof. Let us examine separately the cases that d ti > d t > , d ti < d t > , or d ti — d t ' ■ If d ti > d t >. Si is selected 
irrespective of the trend (left clause). This condition is sufficient because segments cannot cross. If d ti < d t > the 
formula evaluates to false, leading to select S' it which is correct for the same reason. If d ti = d t > then three cases 
are again possible: either both segments are flat (and therefore equals, Si is arbitrarily selected), or both segments 
are a slope (and therefore equals, S[ is arbitrarily selected), or Si is flat and S[ is a slope (a flat segment starting at 
the same high as a decreasing slope will necessarily remain higher; Si is therefore selected), or Si is a slope and S[ 
is flat (symmetrical case, <S- will be selected). □ 

In reality, the segments do cover intervals of various sizes and temporal location, and in some cases the corre- 
sponding distance values may "cross" even when the intervals covered by two segments are aligned (e.g. a slope 
segment starting at a slightly higher value than a flat segment). Rather than complicating the above logic of aggre- 
gation, we pre-process the two tables before aggregation as follows: 

1) Split the segments so as to align both tables: Splitting a segment (ti,dt it trend) at date t\ comes to insert a 
subsequent entry trend) such that d t '. — d ti if trend = flat and d t '. = d ti — (t^ — U) if trend — slope. 
Each table undergoes such a split relative to every index date (i-s) that exists only in the other table. An 
additional split relatively to time may also be added for convenience (if not already present). 

2) Split further to eliminate crossing values: Given two aligned segments of size I, if one is flat and the other a 
slope, say (t, df t ,flat) and (t, d St , slope), and < d Sf — df t < I, then both segments are split at t + (d Sf —df f ) 
precisely. 

Theorem 4. These pre-processing steps produce tables whose segments are aligned and non-crossing. 

Proof. First of all, the fact that a split preserves the consistency of a table is clear from the formulas in step 1 
(i.e., using the same distance value if the segment is flat; decreasing it otherwise by an amount equal to the time 
elapsed since the beginning of the segment). Since each table undergoes a split with respect to every index date that 
exists only in the other table, the resulting tables must contain the same index dates (aligned segments). As for the 
crossing segments, let us first observe that two flat segments cannot cross since their value is a constant, neither can 
two slope segments because their value decrease at a same rate (the rate of time). Therefore, crossings may only 
occur between segments of different trends. Let S = d St — df t be the difference of initial value between the slope 
segment and the flat segment. The segments cannot cross if the slope segment is already below the flat segment at 



the beginning of the interval (a slope is always decreasing), neither can they if S is larger than I because the final 
value of the slope segment will still be above that of the flat segment. Therefore, crossings only occur if < 5 < I. 
Besides, both types of segment trend being linear, two given segments cannot cross more than once. It is therefore 
sufficient to split them with respect to their crossing point t + S. □ 

Figure |7] illustrates the aggregation of two distance tables that relate again to the same example, i.e., d a , t (b) and 
d a ,t (c) in the triangle TVG of FigureH] Note that no situation of crossing segments occurred to be handled. Based on 
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Fig. 7. Aggregation of two distance tables (corresp. to d a ,t(b) and d a ,t(c) in Fig.|4j. 

the final aggregation, the emitter can decide when to initiate the intended fastest broadcast - here, anytime between 
dates 20 (inclusive) and date 29 (exclusive) modulo p. 

5 Concluding Remarks and Open Problems 

The problem of computing fastest journeys in delay-tolerant networks was defined and solved in |4) in its combi- 
natorial variant (i.e., by a centralized algorithm that has full knowledge of the network schedule). We formulated 
this problem in a distributed setting, and solved it in the general case that contacts have arbitrary durations and can 
possibly overlap with each other. This case is more complex than the case with punctual contacts, because both 
direct and indirect journeys can co-exist in the network. Using the T-CLOCKS abstraction from Q, we showed how 
nodes could learn their own temporal eccentricity in such a context in the particular case of periodically-varying 
graph. Based on this information, the task of building fastest broadcast trees reduces to that of building a foremost 
broadcast tree relative to (any of) the date of minimum temporal eccentricity. 

The purpose of this paper was primarily to demonstrate the feasibility of this problem. As far as complexity is 
concerned, most of the communication cost is encapsulated at the level of T-CLOCKS. Characterizing and improv- 
ing this cost thus requires to characterize and improve that of T-CLOCKS. However the complexity of T-CLOCKS 
is still unknown. This open problem is regarded as a relevant research avenue, all the more that T-CLOCKS is used 
as a building block to solve various concrete problems. Contrary to most distributed algorithms for static networks, 
the complexity of an algorithm in the context of interest does not only depend on the number of nodes and edges, 
but is strongly dependent on the number of topological events during the execution (in fact, a vast majority of 
communications and computations are precisely triggered by these events). It seems therefore reasonable to start 
by looking at the way various network schedules could impact complexity. 
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